Strong-field ionization of atoms and molecules: The two-term saddle point method 
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We derive an analytical formula for the ionization rate of neutral atoms and molecules in a 
strong monochromatic field. Our model is based on the strong-field approximation with transition 
amplitudes calculated by an extended saddle point method. We show that the present two-term 
saddle point method reproduces even complicated structures in angular resolved photo electron 
spectra. 
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I. INTRODUCTION 

In order to describe fully the dynamics of molecules 
and atoms subject to an external laser field, one must in 
principle solve the time dependent Schrodinger equation 
including all degrees of freedom. Such ab initio solutions 
are, however, impossible for any but the most simple 
systems, and additionally these methods arc often only 
available for a few specialized theoretical research groups. 
Fortunately much physical insight can be achieved by 
simpler models. For example, many strong-field phe- 
nomena can be successfully interpreted if one uses the 
Ammosov-Delone-Krainov (ADK) tunneling model 1] to 
describe ionization. The success and the analytical sim- 
plicity makes the ADK model ideal for widespread use 
not only for atoms but also for diatomic and poly- 
atomic molecules. 

Along with the ADK model, the strong-field approxi- 
mation (SFA) is one of the most widely used models to 
describe detachment of anions and ionization of atoms 
in intense laser fields. Compared with the ADK model, 
the SFA is more suited for obtaining angular and energy 
resolved spectra. The two models are in fact connected 
since the tunneling rate can be obtained from the SFA 
in the low frequency limit 0, 1^. The initial work by 
Keldysh concerned ionization of hydrogen. The model 
was further developed by Faisal and Reiss 0, 1^ and 
is commonly known as the Keldysh-Faisal-Reiss (KFR) 
model. Later on, the model was extended in various ways 
(see Ref. (2| for a recent review), e.g., to take into account 
rescatterin g long range Coulomb potential in the 

final state Il2j |. multiple electrons [T^ll^ and molecular 
structure Il5| . More systematically, for short-ranged fi- 
nal state interactions the SFA transition amplitude is the 
leading term in an exact ^'-matrix series • 

In order to evaluate the matrix elements that enter 
the expression for the ionization rate in the SFA, one 
may use the saddle point method to obtain approximate 
closed analytical formulas. The saddle point method can 
be applied in both the length- and velocity gauge. In 
the velocity gauge, the saddle point method breaks down 
at intensities below 10^'^W/cm^ Despite its wide- 

spread and long-term use - already Keldysh applied the 
saddle point method in the initial work concerned with 
the limit of small momenta of the outgoing electron and 



Ref. ^3 reviews other limiting formulas - we are not 
aware of a similar study of the limitation of the saddle- 
point method in the length gauge. The main purpose 
of the present paper is to provide a detailed discussion 
of the applicability of the saddle point method in neu- 
tral atoms and to extend the theory to cover molecules. 
In this effort, we identify a straightforward extension of 
the conventional saddle point method. We call the ex- 
tended theory the 'two-term' saddle point method, and 
we show that the present method increases the accuracy 
considerably. 

The paper is organized as follows. In Sec.^we outline 
the theory. In Sec. lIIII we test the accuracy of the saddle 
point method by presenting result on various atoms and 
molecules. Section Hvl concludes. 



II. THEORY 

The saddle point method gives very accurate results 
for detachment rates of negative ions 0| and a sad- 
dle point formula that covers also neutral atoms, irre- 
spectively of the value of the momentum of the outgoing 
electron, is known The application of the latter for- 
mula, however, was not considered until recently ^19.. JCl . 
The theory outlined here follows closely the derivation of 
Ref. . The differences are that (i) we take into account 
molecular structure, and (ii) the previous theory only in- 
cluded one term in saddle-point evaluation of a particular 
integral, whereas we keep two terms to increase the ac- 
curacy and range of applicability [see Eq. H14|) below]. 
Equations IJl-Q) below summarize the basic formulas 
from Ref. 5] and are included here for completeness. 

In the single-active-electron approximation, we con- 
sider the direct transition of the electron in an initially 
bound state to a continuum state \E'q due to the lin- 
early polarized laser field F{t) ~ Fq cos(u;t) with the pe- 
riod T = 27r/w. We quote the expression for the angular 
differential ionization rate (atomic units h = vrie = \&\ — 
1 with the electron charge e = — 1 are used throughout) 
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with the transition amplitude for the n-photon process 



Agn = ^J^ (*g(r,i)|F(t) • r\^o{r,t))dt, 



(2) 



which is to be calculated at the momentum qn = 
^2{nLij — Eb — Up), with Ei, the binding energy of the 
initial bound electron and Up = Fq/{4u!'^) the pondero- 
motive potential. Since the final momentum is real, a 
minimum number of photons n^i^ must be absorbed. In 
the SFA, the interaction between the field and the elec- 
tron in the initial state '^o{r, t) is neglected and accord- 
ingly ^'o(r, t) = <I>o('") exp {lEbt), where $o('") is the sta- 
tionary solution of the field-free Schrodinger equation. 
Additionally, interactions between the residual ion and 
the continuum electron are neglected in the final state 
which is then described by a Volkov wave 



*g(r, t) = exp |z[q + A{t)] ■ r - j [q + A 

(3) 

with the vector potential A{t) = — Fq / lo sa\{uLit) . We 
omit the lower integration limit corresponding to an adi- 
abatical turn-on of the field at i — *■ — oo. 

Following Ref. H we write the transition amplitude 
equivalently as 



where F{a, b; c; z) is Gauss' hypergeomtric series. We 
insert the Fourier transform above in Eq. I^J and obtain 
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X e'^^'Ut, (7) 
with kinematical momentum Q{t) = q + A{t), and action 



S{t) = 



nujt 
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-dt' 



(8) 



, cos{ijjt) — ^ sm{2ujt). 



In a multiphoton process n ^ 1 and the exponential 
factor e*'^'^*) therefore oscillates rapidly on the interval 
< t < T. This fact makes the integral difficult to 
evaluate directly from Eq. (71). 

The time integral of Eq. ([TJ) follows obviously the real 
t axis. From Eq. O and the convergence of the hyperge- 
omtric function, we see that the integrand is an analytical 
function of t except at complex instants of time satisfying 
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with K = \/2E}, and ^o{q) being the Fourier transform 
of ^oir)- The transition into the continuum takes place 
at large distances from the ionic core, and hence, for the 
initial state, it is accurate to use the asymptotic Coulomb 
form which we expand in partial waves 



$o(r) 



(5) 



with V — Zj K and Z the charge of the residual ion. Here 
we assume a general non-spherically symmetric poten- 
tial. For atoms, only one term contributes and for lin- 
ear molecules in the body fixed frame, we only include 
m states corresponding to the projection of the angular 
momentum on the internuclear axis. We determine the 
asymptotic expansion coefficients Cim by matching the 
Hartree-Fock orbital of the most loosely bound electron 
to the form of Eq. lO HJ. The Fourier transform of 
Eq. (HI is 
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X F('l^,^±l;/ + |-5)y..(q), (6) 



[q + A{t)f ^ Qitf 



(9) 



where $o[9 + A{t)] is singular. These points coincide 
with the saddle points S'{t) = in the factor e''^^*^ 
We discuss the complex momentum in appendix ^ and 
refer to Refs. |^ 0| for the analytical properties and 
evaluation of the spherical harmonics on a complex vec- 
tor. The continuation to the complex t plane is straight- 
forward as long as we remember to treat the singulari- 
ties with care. Note that the singularities vanish when 
v = Z — corresponding to detachment of negative 
ions. In Fig. ^ we show the integrand of Eq. (7J) for 
a typical set of laser parameters applied to the ground 
state of hydrogen {I = m = 0, Cqo = 2). The inte- 
gral along the closed contour shown in Fig. ^ is zero 
according to Cauchy's theorem since we carefully avoid 
to enclose the singularities marked by crosses. The inte- 
grand is clearly invariant under the periodic translation 
Re{t) + ilm{t) Re(< -I- T) -I- ilm{t). Hence, the contri- 
butions to the integral from the vertical paths V2 and V4 
cancel exactly and, consequently, the contributions along 
the horizontal paths must also cancel. We can therefore 
equally well evaluate the integral Eq. (TJ along the nega- 
tive V3 path. At the expense of introducing complex and 
unphysical times, we see from Fig.nthat we can calculate 
the integral more efficiently in the complex t plane. It is 
apparently much easier to evaluate the integral along the 
path —Vs than along the real axis Pi. Along the former 
path, the factor e*'^'*^ is nearly zero everywhere except 
at the two saddle points while the same factor oscillates 
along the entire path Vi . This fact is of course the basis 
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FIG. 1: (Color online) Integration contour. Upper (lower) 
panel: real (imaginary) part of the oscillating factor e''^'*\ In 
the regions bounded by dashed squares, function values are 
multiplied by 2 x 10^ {Im{Ljt) > 1.5). The parameters corre- 
spond to ionization of hydrogen in the polarization direction 
by 10 photons at 800 nm at an intensity of 1 x 10^^ W/cm^. 



for any asymptotic expansion. Strictly speaking, the sad- 
dle point method requires that we deformate the path to 
pass across the saddle point in the direction of the steep- 
est descent. In practice, the required deformation from 
the horizontal line is small and has negligible influence 
on the final results. 

In the saddle point method outlined in Ref. 0, one 
neglects the variations of Q{t) over the range where the 
factor [Q{ty+K'^]~'^e^^^*'> has a significant amplitude, i.e., 
according to Eq. we let Q{t) = ±iK in the remaining 
factors and obtain 

\n « -5lG'™r(z/ + i) (-)'' ^ (±i)'rz„(q^) 
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where (j) = tot and S'{(j)) = [Q{t)'^ + k'^]/{2uj). We expect 
such an approximation to be most accurate for I = since 
the factors QityYim[Q{t)] in Eq. are constant in this 
case. The integrals are to be evaluated along the negative 
direction of the path V3 near the /^'th point of stationary 



phase. In Ref. |^, the denominator was expanded as 

[^'(</-)]-''«r(0^)]-''(0-</.,,)-^ (11) 



near the saddle points. As we show in Sec. IIIlAI below, 
we obtain higher accuracy if we expand S"(0) to second 
order around the saddle point 
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and by the first order binomial series 
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With this expansion inserted in Eq. (|10|l . the integral now 
becomes a sum of two terms 



with the conventional saddle-point term 

giS(0) 
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and the present correction term 
^ Xo 2[5"(0^)]-+i 
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In Eqs. (|15|l and 1)16(1 . we extended the integration limits 
to infinity and used the asymptotic approximation 
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Jqh recovers the result of Ref. (51 while is a correction 
that arizes from the higher order expansion of [5"(0)]~'^. 
We present the main formula of the present work in the 
next equation 



Im A'=l,2 ^ 
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With the inclusion of Io^{l +Cp), we refer to Eq. (|18|l as the two-term saddle point approximation while neglecting 
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Cf^ and maintaining only Xo^ in l|18() is referred to as 
the one-term saddle point approximation. We note from 
Eq. 1)16(1 that there is no significant additional numerical 
complications involved with the inclusion of the second 
term. 

We find <j>fi from the saddle point conditions of Eq. 



S' = n - z + Csin0^ + 2zsm^ = 0, (19) 

with ^ = —Fo ■ q/uj'^ and z — Up/ui. The solutions to 
Eq. ^ are 
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from which it follows that 
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±^Jl-siTi^(t>^, (21) 

(cos 4>f, + i sin 0^ )"e-* '^^(«+^ "^"(22) 
cos(/)p(^ + 4zsin0^) (23) 



= — ^ sin ( 



42 [2 cos^ 
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where the signs correspond to Q(<^) = T^'^ &t the saddle 
points. We combine Eqs. and it^ - lt^ 

to obtain the analytical approximation to the transition 
amplitude, Eq. (|7|). 

The inclusion of nuclear motion in the molecular case 
was discussed in detail for ionization pH and for high 
harmonic generation within the SFA. The form of 
the amplitude (0 stays the same and the formulas are 
straightforwardly generalized. When it comes to the 
assessment of the accuracy of the saddle-point method 
which is the main objective of the present paper, nuclear 
motion is unimportant and is left out for clarity. 



III. RESULTS 
A. Test case: atomic hydrogen 

First we consider ionization of ground state hydrogen. 
We use this system to benchmark the accuracy of the 
saddle point method against numerical integration. The 
atomic structure parameters are Coo = 2, k = 1 and 
1/ = 1. Furthermore, the asymptotic form of Eq. |SJl is 
identical to the exact wave function at all distances and 
the Fourier transform is 

The spherical harmonic and the hypergeomtric function 
in Eq. {Tj) are both constant and it is therefore exact to 
neglect variations therein around the saddle point when 
we derive Eq. (fTH)! . 

After choosing the alternative integration path —Vz of 
Fig. n the transition amplitude reduces to 
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FIG. 2: (a) Real and (b) imaginary parts of the integrand 
giS(0)|^/(-^-)j-^ along the integration path -Vs in the neigh- 
borhood of the second saddle-point in Fig. solid. The long 
dashed line is the approximation around the second saddle 
point including only Xqm- The two-term approximation in- 
cluding Xop(l + C^) is the short dashed line which overlaps 
completely the solid line. The parameters of the laser are as 
in Fig. [U 



Here we test the accuracy of Eq. (fTH|) [or and in 

particular the difference between the one-term (only Xo/i 
included) and the two-term (Jo/i(l+C^) included) saddle 
point formulas. In Fig.|21we show the integrand along the 
integration path in the neighborhood of the second sad- 
dle point (see Fig^. Additionally, we show the results 
of the approximations using the one- and two-term ex- 
pansion around the second saddle point in Eq. ((13(1 . We 
see that the one-term expansion recovers quite well the 
peak structure around the saddle point. In the wings of 
the peak, the two-term expansion is significantly better, 
which is most evident from the real part of the integrand, 
panel (a). We have integrated the integral numerically 
and obtained the value 5.64 x 10~^ while we obtain the 
values 4.55 x 10^^ and 5.39 x 10^* for the one- and two- 
term saddle point approximation, Eq. 1(14(1 summed over 
the two saddle points. 

Having seen that the saddle point method is accurate 
in the single case above, namely ionization parallel to the 
field by 10 photons, we now present the n-photon angular 
differential rates at varying photon orders in Fig.O Here 
9q is the polar angle of the outgoing electron with respect 
to the polarization axis. Figure |2| shows that both sad- 
dle point methods predict an angular structure in close 
agreement with the numerical integration. The rates ob- 
tained by the single-term approximation are, however, 
around 35% too small. The two-term approximation is 
significantly better with an accuracy within 10%. 

In our final test, we consider the total ionization rate 
integrated over all angles of the outgoing electron. In 
Fig- El (a) we present total ionization rates at 800 nm ob- 
tained with the numerical integration and the one- and 
two-term saddle point approximation. All three methods 
produce results in quite good agreement over many orders 
of magnitude on the scale shown in the figure. In order 
to investigate the accuracy of the saddle point method in 
some more detail, we calculate the ratio VFsaddic/Wnum 
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FIG. 3: (a)-(d) Angular differential ionization rate of hydro- 
gen for the lowest number of photon absorptions, n = 10 — 13 
ordered by increasing n. The solid curve is obtained by nu- 
merical integration while long and short dashed curves are 
obtained by the saddle point method with the one- and two- 
term approximation, respectively. The laser wavelength is 
800 nm and the intensity is 1 x 10^^ W/cm^. 
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FIG. 4: (a) Absolute rates at a wavelength of 800 nm. The 
result of numerical integration is indicated by the solid line 
while the long- and short dashed lines are the one- and two- 
term saddle point approximation, respectively, (b) and (c) 
Ratio between total ionization rates obtained by the saddle 
point method and numerical integration for varying wave- 
length and intensity. We show saddle point results for the 
one- and two-term approximation in panels (b) and (c), re- 
spectively. 



Ionization of atoms 



between the rates obtained by the saddle point method 
and the numerical integration. Figures 01 (b) and (c) 
present the results for various wavelengths and intensi- 
ties. Again, we use both the one- and two-term approxi- 
mation, i.e., we study the results of Eq. (|18|l withZo/^ ^nd 
2^0At(l respectively. First, we see that the accura- 

cies of both saddle point methods are nearly independent 
of the intensity for each fixed value of the wavelength. 
From panel (b) we note that the results of the one-term 
approximation depend significantly on the wavelength. 
The error is up to a factor of two for the shortest wave- 
length 248 nm while the error decreases with increasing 
wavelength. The two-term approximation, on the other 
hand, produces much more accurate results, panel (c). 
The rates are approximately 10% too small for all in- 
tensities and wavelengths considered. Even though the 
simple single-term saddle point approximation is some- 
what poorer than the two-term approximation the error 
in the hydrogenic case is approximately constant over 
twelve orders of magnitude and is not expected to be of 
major significance compared with the approximations in 
the SFA itself. In the rest of the paper, we use only the 
two-term saddle point formula. 



In this section we show results for the noble gas atoms 
where the active electron initially occupies an orbital in 
the filled p shell. We calculate the rates for each of the 
states m = — 1, 0, 1 and multiply the result by two corre- 
sponding to two equivalent electrons in each orbital. We 
take the atomic structure parameters Cim and Ef, from 
Ref. 

In Fig. El we present the absolute and relative rates for 
the argon atom. As in the case of hydrogen, the saddle 
point method is accurate over many orders of magni- 
tude. Fig. 0(a). Interestingly, the saddle point method 
is slightly better for the m = 1 [panel (c)] state compared 
with the m = state [panel (b)]. This m-dependent ac- 
curacy turns out to be important in the molecular case 
as we show in Sec. IIII Cl below. 

Figure shows angular differential rates summed over 
all photon absorptions at a wavelength of 800 nm for dif- 
ferent intensities. We show the results for both the m = 
and m ~ 1 state and again we see that the saddle point 
method works better for the m — 1 state. The general 
features in the angular spectra are, however, very well 
reproduced in all cases. 

We mention in closing that the results for krypton and 
xenon are very similar to argon and are therefore omitted 
here for brevity. 
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FIG. 5: (a) Absolute rates for argon at a wavelength of 
800 nm. The results of numerical integration are indicated 
by the solid lines while the short dashed lines are the two- 
term saddle point approximation. The upper and lower sets 
of curves are for m = and m = 1, respectively, (b) Ratio be- 
tween total ionization rates from the m = state obtained by 
the saddle point method and numerical integration for vary- 
ing wavelength and intensity, (c) Similar to (b) for the m = 1 
state. 
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FIG. 6: Left (right) column: Angular differential rate for 
argon in the m = (m = 1) state. Solid lines are obtained by 
numerical integration while the dashed lines are by the two- 
term saddle point method. The laser wavelength is 800 nm 
and the intensities are (a)-(b) 1 x lO^^W/cm^, (c)-(d) 1 x 
10''' W/cm^ and (e)-(f) 3 x lO'"* W/cm^ 



C. Ionization of molecules 

In the molecular case, the calculations are most conve- 
niently performed in the laboratory fixed frame with the 
z axis parallel to the laser polarization. Accordingly, we 
must express the initial wave function in this frame. The 
wave function and asymptotic expansion coefficients Cim 
are, however, most naturally expressed in the body-fixed 
molecular frame. If the body-fixed frame is rotated by the 
Euler angles (a, /?, 7) with respect to the laser polariza- 
tion, we rotate the wave function into the laboratory fixed 
frame by the rotation operator ^oif) D{a, (3, 7)<I>o(r). 
The rotation operation effectively allows us to express 
the asymptotic coefficients in the laboratory frame (LF) 
by the corresponding coefficients in the molecular frame 
(MF) 



I 

E 



(27) 



where 



^(0 



(a,/?, 7) is a Wigner rotation function 



|2J| . For linear polarization and the linear molecules con- 
sidered in the present work, we only need to consider 
rotation by the angle /3 between the molecular and field 
axes. We refer to Ref. for the coefficients Cf^ . 

Figure H presents angular differential rates for differ- 
ently aligned N2 molecules which ionize from the dou- 
bly occupied Scr^ orbital. We show both numerical [(a)- 
(c)] and two-term saddle point results [(d)- (f)]. First we 
note that the two methods agree perfectly on the shape 
of the angular distribution for all alignment angles, 0. 
The structures are also in good agreement with Ref. [23 , 
where we used an atomic centred Gaussian basis expan- 
sion for the initial state and calculated the transition 
amplitude numerically. Secondly, the overall structure is 
nearly independent of /3. The angular rate is simply much 
favoured along the polarization direction in all geome- 
tries. This observation agrees well with the predictions 
of tunneling theory where the electron by assumption es- 
capes along the polarization axis (For Kcldysh parame- 
ter 7 = KUJ / Fq <C 1 , the ionization dynamics is tunneling 
like. In Fig.[7| 7 — 0.81, i.e., approaching the tunneling 
regime.) 

Figure |H1 shows the angular differential rate for O2 
which ionizes from the two half-filled degenerate tt^ or- 
bitals, one with m = 1 and one with m — —1. We show 
the results for a single electron with projection m — 1 and 
note that the rate is similar for m = —1. Again we see 
that the two methods predict the exact same complex an- 
gular structures. The structures can be understood from 
the symmetry of the initial wave function. The initial TTg 
orbital has zero amplitude along- and in the plane per- 
pendicular to the molecular axis and the nodal structure 
of the wave function forbids the electron to be emitted 
along the vertical polarization axis when this axis coin- 
cides with a nodal plane 3]. 

In Figs. [7|and|Slwe see that even though the angular 
structures agree perfectly, the absolute scales differ by up 
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(a) (b) (c) 




FIG. 7: (Color online) Angular differential ionization rate for 
N2 aligned at an angle of 0° [(a) and (d)], 45° [(b) and (e)], 
and 90° [(c) and (f)] with respect to the polarization. Panels 
(a)-(c) are obtained by numerical integration and panels (d)- 
(f) by the two-term saddle point method. In all panels, the 
polarization direction is vertical. In panels (a) and (d) the 
molecular axis is along the vertical polarization. In the other 
panels the orientation of the molecular axis is indicated by the 
line through the origin. All three coordinate axes are scaled 
equally and the rates are given in units of s~^ on the scale 
indicated in each panel. The laser wavelength is 800 nm and 
the intensity 2 x 10^* W/cm^. 
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FIG. 8; (Color online) Similar to Fig.[7|but for O2. 



to a factor of four. We therefore finally turn to a discus- 
sion of the alignment dependent rate. We calculate the 
total rate integrated over all angles of the outgoing elec- 
tron and show the results in Fig.|51 The figure shows that 
both methods agree that the rate for N2 is maximized 
when the molecule is aligned parallel to the polarization 




45 90 135 180 45 90 135 180 

/3 (degrees) (degrees) 

FIG. 9: Alignment dependent ionization rates for (a) N2 and 
(b) O2. The solid lines are obtained by numerical integration 
and the dashed line by the two-term saddle point method. 
The parameters of the laser are as in Fig. 13 



axis (/3 = 0) and minimized when aligned perpendicu- 
larly {f3 — 90°). Such an ahgnment dependence is also 
seen experimentally . For O2 both methods also agree 
that the rate is maximized around an alignment angle of 
40°. For both molecules, however, the two-term saddle 
point method predicts a too large variation compared 
with the numerical integration. The reason for this dis- 
agreement lies in the m-dependent accuracy of the saddle 
point method as we discussed in Sec. IIIIBI When we ro- 
tate a molecule we change the expansion of the initial 
wave function in the laboratory fixed spherical harmon- 
ics. The rotation operation mixes the different m-states 
according to Eq. (|?7jl . Since the part of the transition 
amplitude that corresponds to each of the m-states can 
be either slightly too large or too small by the saddle 
point method, the overall accuracy depends on the partial 
wave decomposition after the rotation. For the atomic 
I = 1 and m — 0, ±1 states of Sec IIII Bl the two-term 
saddle point method is still quite accurate and the small 
differences reported there compared with the numerical 
integration cannot account for the disagreement between 
the two methods in Fig. For the molecules, however, 
we have included angular momentum states up to Z = 4 
and it turns out that the saddle point method becomes 
increasingly inaccurate with increasing I. If the active 
electron initially occupies an orbital with a component of 
non-zero angular momentum, we expect the saddle point 
method to be somewhat poorer than for Z = as we dis- 
cussed in deriving Eq. (|10|l . In Eq. (|10|l . we evaluate the 
factors Q{ty Yijn[Q(t)] and the hypergeomtric function at 
the saddle points. This approximation is naturally most 
accurate if Q{t) is nearly constant in the vicinity of the 
saddle points. In appendixIXl we calculate Q{t) along the 
integration contour and we see from Fig. llOl that the vari- 
ation in Im[Q(i)] is in fact close to maximal at the saddle 
points. If we require higher accuracy of the saddle point 
method, we must take at least the first order variation in 
Q{t) into account and modify Eq. H10|) accordingly. 



8 



IV. SUMMARY AND CONCLUSION 

Based on the length gauge SFA, we proposed a two- 
term saddle point formula which is applicable to neutral 
atoms and molecules. We presented calculations on var- 
ious atoms and molecules with the primary aim to test 
the accuracy of the method. The two-term saddle-point 
evaluation is very accurate in the case of ionization of hy- 
drogen while the accuracy is within 10% for noble gas 
atoms which undergo ionization from a p shell. Remark- 
ably, the structures of the angular photo electron spec- 
trum predicted for all systems are in perfect agreement 
with numerical calculations. We have identified that the 
saddle point method in our formulation works best if the 
initial wave function is a zero angular momentum state. 
Multicentric molecular wave functions contain many dif- 
ferent angular momenta and correspondingly we see small 
inaccuracies when we use the saddle point method for 
molecules. 

In contrast to previous reports on saddle point meth- 
ods in the velocity gauge SFA , we did not find a criti- 
cal lowest intensity below which the saddle point method 
fails. Even though we find small errors which are direct 
consequences of using the saddle point method instead of 
a numerical evaluation of the transition amplitude, the 
error is nearly constant for a wide range of intensities and 
is small compared to the large variations in the absolute 
rates. Furthermore, we should keep in mind that the SFA 
itself is only the leading order of an S'-matrix series. The 
small error in the saddle point evaluation may therefore 
turn out to be insignificant compared to, e.e^ neglecting 
the Coulomb interaction in the final state |27l |. 

We conclude that the saddle point method in the 
present two-term version can be used with advantage for 
long wavelengths and high intensities when many photon 
absorptions lead to ionization. In this case, the transi- 
tion amplitude is difficult to evaluate numerically since 
the integrand oscillates rapidly. The SFA also applies to 
non-monochromatic fields, e.g., a few cycle pulse. The 
transition amplitude is then calculated by an integral 
over the full duration of the pulse. Numerical integration 
by standard Gaussian quadrature requires thousands of 
function evaluations to obtain convergence [2^ . It would 
clearly be desirable to extend the saddle point method to 
such a situation where we need just a few saddle point 
evaluations. 




(b) 









-1 n 1 



7r/2 IT 37r/2 27r 
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FIG. 10: (a) Q{t)^ in the complex plane along the path —P3 
in Fig. (b) Real (solid) and imaginary (dashed) part of 
Q(t) according to the definition Eq. I|A3|I . 



APPENDIX A: COMPLEX MOMENTUM 



In connection with Eq. 
matical momentum 



Qit) 



lO, we introduced the kine- 



Fo 



siii{ujt). 



(Al) 



We let the laser polarization point in the z direction and 
find the squared momentum for a complex time t = tr + 
iti 



Qitf = ql + ql + [qz 



— sin(t^t)]2 (A2) 

UJ 



with sin^Lut) — siii{u!tr) cosh{ujti) + i cos{u!tr) smh{Luti) . 
We wish to calculate Q{t) along the path in Fig. ^ 
The path is parametrized as {t — tr + iti\0 < tr < 
T} with the imaginary part ti fixed. In the polar form 
Q(t)2 = \Q{t)'^\e^^^*\ we define the domain of the phase 
of Q{t)^ between — tt < < tt. When we calculate the 
square root, we lie a branch cut along the negative semi- 
axis and change the sign when the branch cut is crossed. 
In Fig. ^l(a), wc show a parametric plot of along 
the path —V3 from Fig. ^ Using the definition above, 
we find Q{t) 



Qit) 



+ v /|Q(^)2| e»^W/2 Outer loop 
-^|Q(t)2|e*^(*)/2 Inner loop 



(A3) 
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where the outer and inner loops refer to Fig. EH (a). In 
Fig. ^1 (b), we show the real and imaginary parts of 
Q{t) along —V3. As in Fig. ^ we have k = 1 for the 
ground state of hydrogen. We see that Q{t) = —in at 
the left saddle point {ujtr « 0.29) and Q{t) = in at the 
right saddle point {uitr ~ 2.85), which leads to the factor 
(±1)' in Eq. (P|l . 
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